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ABSTRACT 

Research  on  wind  and  temperature  profiles  over  land  have 
established  a  number  of  fundamental  relationships  that  are  tested 
here  over  an  oceanic  location.     Most  attempts  at  calculating  the 
roughness  parameter  over  oceanic  surfaces  have  been  based  on 
surface-layer  theories,  and  often  employ  an  iterative  procedure  in 
approaching  the  sea  surface.    This  paper  employs  a  similarity 
relationship  between  the  wind  and  temperature  profiles,  using  data 
at  ship  "Victor"  appropriate  to  forced  convective  theory,  as 
defined  by  Priestley. 

The  theory  of  forced  convection,  as  applied  here,  makes  it 
necessary  to  compute  the  Richardson  number  and  the  Monin-Obukhov 
scale  length  as  preliminaries  to  the  computation  of  roughness 
parameter.    The  roughness  parameter  and  friction  velocity  are 
obtained  from  the  theory  using  an  iterative  procedure.    The  wind 
speed  at  20  meters  is  correlated  against  scale  length  and  roughness 
parameter.    The  latter  makes  very  little  contribution  to  the  variance 
of  the  windspeed  while  the  scale  length  contributes  in  a  manner 
reasonably  expected  from  theory. 
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List  of  Symbols 

Symbol  Definition 

u  Wind  speed  at  n  meters  height 

Z  Height  in  centimeters  above  the 

water  surface 

k  Von  Karman    Constant 

u  Friction  velocity 

L  Monin-Obukhov  scale  length 

* 
T  Scaling  temperature 

H  Eddy  jhjfeat  transport  in  the 

vertical 

Rj  Richardson  number 

g  Acceleration  of  gravity  (980  cms) 

Cp  Specific  heat  of  air 

P  Density  of  air 

Z0  Roughness  parameter 

Q  Potential  temperature 

S  Non-dimensional  wind  shear  after 

Monin  &  Obukhov 

Kjj  Eddy  diffusivity  for  heat  conduction 

K»,  Eddy  diffusivity  for  momentum  flux 


Kg_2Q  Ratio  of  wind  at  600  cms  to  that 

at  2,000  cms. 


iv 


Symbol  Definition 

K     .  Ratio  of  (Z |^  )j  4  to  the  wind 


1.4 


difference  at  20  meters  minus  that 
at  ten  centimeters 


"Z>  Surface  stress 

Cj   4  A  more  exact  form  of  K^  4 

Vh/fZQ  Surface  Rossby  number 


1 .    Introduction 

Studies  of  the  roughness  parameter,  ZQ(    have  been  numerous 
over  land  and  a  number  of  fundamental  relationships  using  various 
meteorological  parameters  have  been  determined.    Far  less 
research  has  been  accomplished  in  this  field  over  the  oceans. 
This  is  readily  understandable  when  one  considers  the  expense  and 
difficulty  in  erecting  and  maintaining  an  instrumented    I  platform 
at  sea.    With  this  in  mind,  one  could  turn  to  what  should  be  a 
prime  source  of  data — the  ocean  station  weather  ships. 

Certain  problems  become  evident  very  quickly  when 
considering  station  Victor.    Four  different  ships  occupy  this 
station  and  only  two  have  the  same  anemometer  height.    This,  at 
first,  seemed  to  present  an  insurmountable  problem  in  using  the 
data  but  the  theoretical  error  involved  in  using  a  median  height  of 
twenty  meters  instead  of  actual  heights  is  less  than  two  and  one 
half  percent.    When  considering  that  winds  are  recorded  only  to 
the  nearest  knot  and  temperature  to  the  nearest  degree,  the  error 
introduced  in  the  wind  profiles  by  considering  all  heights  equal  to 
twenty  meters  is  insignificant. 

Monin  and  Obukhov    f3  J  /  in  their  universal  similarity  theory 
introduced  three  scaling  parameters  that  are  considered  essentially 
invariant  with  height  for  layers  near  the  ground.    These  three 


parameters  are  scaling  velocity,  u*„  scaling  length,  L,  and 
scaling  temperature,  T  ,  and  have  been  well  described  in  a 
recent  text  (Lumley  and  Panofsky)     flj  on  atmospheric  turbulence. 


L,  u*  and  T*  are  defined  as  follows: 


L  =  -  tfffi1'  (i) 


1     =  -  J?u*cPP  (2) 


U*  -     \I—Z —  (3) 


It  should  be  noted  that  a  selection  of  cases  has  been  limited 
to  cases  for  which  -RS^.  032  holds.    According  to  Priestley's  dri^i- 
rion  DfCJ   our  restriction  limits  the  study  to  the  cases  of  "forced 
convection".     Here  Ri  refers  to  the  turbulent  parameter 


(4) 


and  is  to  be  distinguished  from  the  flux  Richardson  number  Rf^ 
defined  as 


(5) 


Under  normal  circumstances,  it  is  considerably  easier  to  deal 
with  Ri  than  with  Rf ,  since  the  factor  KH  /  KM  is  itself  an 
increasing  function  of  Ri. 
2.      Important  properties  of  the  forced -con  vective  surface  layer 

According  to  the  similarity  theory  for  the  near-neutral,  or 
forced  convection  case  (as  defined  in  section  1),  the  wind 
profile  can  be  written 

In  a  like  manner  the  temperature  profile  can  be  written,  after 
Lumley  and  Panqjfofey 


*-a=tTa.£+{?) 


(7) 


In  both  (6)  and  (7),  Lumley -andrPanofafey  CJ   suggest  the  value 
P'  =4.5. 

By  dividing  equation  (7)  by  equation  (6)  and  taking  the 
derivative  of  both  numerator  and  denominator  with  respect  to 
height,  one  has  a  relationship  governing  the  constants  T*  and 
u*/k: 


ign  -  m/b-z-  (8) 


Relationship  (8)  means  that  the  change  in  the  temperature 
over  a  given  height  is  proportional  to  the  change  in  the  wind  over 
the  same  height.    Even  more  important,  this  ratio  is  constant 

*  * 

with  height  since  T   ,  u   ,  and  k  are  all  considered  constant  in  the 
surface  layer. 


The  ratio  of  these  constants  T*/ (u*/k)     may  be  determined 
from  the  air  temperature  at  psychrometer  height  (six  meters)  and 
the  temperature  at  ZQ  ,  assumed  to  be  the  same  as  the  sea 
surface  temperature,  so  that  equation  (8)  may  be  written,  in 
finite-difference  form  centered  at  three  meters  as 


fit)    ~"    Ui~U° 


(9) 


Here  uQ  is  the  wind  speed  at  Z0  and  is  zero.  By  equation  (8)  this 
ratio  is  constant  at  all  levels  in  the  surface  layer,  including  that 
at  1.41  meters. 

From  the  definition  of  L,  it  follows  that  L  satisfies  the 
identity 


e=fefos 


(10) 


However,  Priestley/£4¥has  introduced  a  surface  layer  scale 
lenght  L*  given  by 


L=Lft- 


(ID 


so  that  (10)  becomes 

jf  •  =  R;  S  (12) 

a  formula  lacking  the  troublesome  ratio  Kjj/Km,  which  has  been 
found  to  depend  upon  Ri.     But  from  the  definition  of  T*,  we  have 


or 


is*r-=ltftj 


.Kh  fyp£  -  -**•$&**  <13» 


From  (13),  together  with  the  similarity  profiles  (6)  and  (7),  which 
apply  in  the  case  of  forced  convection,  it  follows  that  Ku=K     , 
and  therefore  L=L' .    This  is  a  fundamental  property  characteristic 
of  similarity  profiles  in  forced  convection. 
3.      Derivation  of  a  specific  value  for  (Ri)i#4i 

From  equations  (4)    and  (9),  the  Richardson  number  is 
expressible  in  the  form 

R.  _  %fc&^  (z^Uz);.//(  <14> 

where  every  parameter  in  (14)  is  now  known  except  for  the 
quantities  u6  and  (  Z§2.  )\.A' 


The  surface-layer  wind  profile  is  assumed  to  be  valid  to 
heights  of  twenty  meters.    The  windspeed  Ug  may  then  be  related 
to  U20  by  a  new  parameter  through 

U4  =  U»  Kws  (15) 

where  Kg_20  is  formed  from  the  ratio  of  Ug  to  ^g.    This  now 
leaves  the  quantity  {2L  flZ    )i    41  as  tne  so^e  unknown  in 
equation  (14) .    One  may  approximate  this  quantity  by  forming  the 
ratio  of  two  wind  differences  centered  at  1.41  meters.    In  the 
following  equations,  subscripts  denote  the  height  of  a  variable 
in  meters  although  substitution  in  centimeters  is  generally 
required.    For  example  we  have 


and 


Ml^j.  in  iS2 


z 

XTToo  M  A'         CJLzoo, 

Therefore,  it  follows  that 


ifcr^~  =  Uzo  K... 


(17) 


(18) 


Now  the  last  two  unknowns  in  equation  (14)  arefoUnd  in  terms  of 
known  or  measured  quantities.    Following  Martin  £2^  , 


1/     -   Lldo^^uj 


and,  similarly, 


(U  ¥  +  i2QQ) 

^o  -  l  -  j  a  L  +  gfesr  (20) 

with  all  values  in  centimeters.    Kg_2n  may  be  approximated 
satisfactorily  for  obtaining  Ri  by  omitting  In  Z0. 
The  right  sides  of  (l5)  and  (16)  are  now  known  as  a  function  of  L 
so  that  (14)  may  be  solved  for  Ri  as 

>W  K,,A.)!v-;    (D    SU£  (21) 


where  K^   4,  and  Kg_90  are  given  as  functions  of  L  in  (19),   (20). 

In  the  forced  convection  regime,  Lumley  and  Pan  of  sky  show 
that  the  Richardson  number  may  be  written  as 

Ri  -  TT   iez        >     KH  -  K^  (22) 

Eliminating  Ri  between  (21)  and  (22)  leads  to  a  quadratic  in  L, 
the  only  meaningful  root  of  which  is  given  by 

L  r    A  +  ^-O* 


A  =  \.,\'\(JLZ0O)-.7ln%-tH)OO  Ri'(JL20O+L20OO) 

B  -   2.  ft  Ri'(Xv,  2 0O)(JU  2 000)  (24) 

c  =  4D:tRi'f(i-20O)(/L  zaxw\*.\d) 

-Wf  Rt'CL  2  0D)(L200O)J 

where 

K'  r  tufT 


Ri  may  now  be  obtained  by  solving  equation  (22)  with  Z=141.1 
cm  and  the  value  of  L  obtained  from  (23). 

The  quadratic  equation  defined  by  (23)  and  (24)  had  as 
input  data  only  windspeeds  at  20  meters,  and  temperature  - 
differences  measured  from  the  surface  to  6  m.    Only  one  data-set 
was  selected  at  ship  Victor  per  day,  and  these  were  selected  with 
a  view  toward  insuring  cases  of  forced  convection:  for  example, 
a  vertical  temperature  difference  less  than  or  equal  to  6C,  and 
windspeeds  as  large  as  possible  under  these  circumstances  were 
selected.    A  table  of  Ri  values  computed  as  a  function  of  both  U20 
and    &&' &o     is  included  in  the  Appendix. 


4.       Iterative  procedure  of  solution  of  roughness  parameter 

Following  Martin C 2  "2,  use  is  made  of  equation  (6)  to  obtain 
the  quantity  U2g  -  Uq    i  in  the  form 

^-U.(  =  ^(2.J02^30l03)+»-2PXV)         M 

But  with  U2q  as  follows 

Uv-^ozUtf)^  +  (OT 

one  obtains  the  result,  U0  j,  given  in  (27) 

U.,  =. 005 U^  +  if (2. 2"? 88 XI -l^&JL^  Zo)  <"> 

By  using  the  Monin-Obukhov  normalized  wind  shear 

O  -   A     /    G>U, 

and  (12),  it  follows  that 

if  =  (2U,  W,  =  flfcflft?*, 

Because  Z0   is  small  compared  to  20  meters,  the  quantity 

(  ^U/^J'viZ  )         will  be  rewritten  involving  a  larger  span,  as 

follows 

Uz-y, 


where  C^   4  will  be  defined  immediately  below. 

(AIM' *  ^  =G,(u*-u*0 


From  (6),  we  the.Q.-have 


Cm  =  Lfe^/(Ute-U.i)-  ^^}'I  uJM  -  (29) 

Combining  (28)  and  (29)  leads  to 


if  =  P#LC.,(LlM-Uw) 


where  C^^  is  defined  in  (29). 

A  third  equation  used  in  the  iteration-process  is 


Uo.,=  ^(£.302^-^,^X1+  J£)  <31> 

but  this  has  already  been  built  into  equation  (25).    The  latter  may 
be  used  for  more  accurate  determination  of  Z0  once  an  earlier 
iterate,  or  satisfactory  guess,  has  been  established.    Equation 
(31)  assumes  the  use  of  equation  (6)r  to  an  elevation  of  10  cnf  above 
mean  sea  level. 

To  adopt    an  efficient  iterative  method    which  does  not  lay- 
too  much  emphasis  on  the  last  10  cm,    we  must    give     proper 
weight  to  the  layer  0-20  m.     Here,  we  use  equations  (27)  and  (30) 
with  u0   1  eliminated.    This  leads  to  the  result 

10 


if- Q,  (V'l  faWJ 


CmT^J,,  (2  2988-2.2171X^7 


(32) 


and,  of  course,  we  also  have 

U2a  r  £  (U  200  +  jU  ^  -f- 1^0 )  03) 

The  iteration  process  may  be  summarized  as  follows: 

1.  Assume  ZQ'  '  =  0.1  and  obtain  a  corresponding  estimate 
of  (u*A)(1)  from  (32). 

2.  With  this  first  estimate  of  (u*A)       ,  use  equation  (33), 
with  U20  and  L  known,  to  determine  ZQV  '. 

3.  Repeat  the  procedure  using  Z0'2'  in  (32)  to  get  a  second 
estimate  of  (u*A)^2^  . 

4.  Repeat  the  iteration  until  a  final     residue   is  obtained: 


When  this  has  been  accomplished  ZQ  is  taken  to  be    Z         ,  and   (tf/k) 

is  the  appropriate  value  from  either  (32)  or  (33). 

One  safeguard  must  be  followed  in  the  use  of  (32):  the 

denominator  of  the  right  side  must  not  be  allowed  to  become  zero 

or  negative.     However,  the  use  of  (33)  helps  to  restrict  (u*/k)  to 

reasonable,  non-negative  values. 

11 


N 


The  results,  giving  ZQ,  u*/k  as  well  as  U2q#  Ri'  /  L  and  Ri 
are  listed  in  Table  1  (Appendix).    In  addition,   some  statistical 
relationships  between  the  parameters  are  also  considered  in  the 
following  section.    The  raw  data  considered  was  drawn  from  30 
randomly  selected  dates.    The  value  of  L  is  computed  from 
equation  (23),  while  Rij.41  is  obtained  from  (22). 

The  computations  were  made  using  data  cards  obtained  from 
the  National  Weather  Records  Center,  Asheville,  North  Carolina, 
for  ocean  station  "Victor".    All  cards  selected  were  for  0300Z,  the 
observation  nearest  noon  local  time,  from  the  months  of  November 
and  December  of  1960  and  of  March  1961.    The  number  of 
intervals  between  successive  days  was  considered  a  random 
function  of  the  date.    Wind  values  were  tabulated  in  this  paper  in 
centimeters  per  second,  and  temperatures  in  degreejCelsius.     Sea 
conditions  during  the  period  were  not  considered, with  the 
exception  that  no  major  storms  were  in  the  immediate  area. 
5  .     Result<s<s3irod3 GQttClU(SiilQfoS.3.lOn 3 

The  values  of  U2Q/  L  and  ZQ  for  the  sample  of  30  cases  were 
subjected  to  a  multiple  regression  analysis  leading  to  a  result  of 
form 

LU    ^-30M&Z«-O.OOS3lLt757el7       (so 

variable  3  variable  1  variable  2 
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Of  the  two  independent  variables  tested  in  (34),  ZQ  plays  an 

insignificant  part  in  that  its  partial  regression  coefficient  was 

-0.00969.    When  the  mean,  value  of  ZQ  =  0.1066     cm  was  inserted 

into  (34),  the  simple  regression  equation 

(35) 

bUo  =  -0.OO53IL  +75H-.74 

resulted.    The  variable  L  had  a  partial  correlation  coefficient  of 
-0.56153,  so  that  its  effect  upon  the  reduction  of  the  variance 
of  u20  was  sufficiently  strong     that    resulting  F-value  due  to  the 
regression  involving  L  was  6.576.    This  F-value,  with  2  and  27 
degrees  of  freedom,  is  significant  at  a  level  in  excess  of  the  99 
percent  confidence  level.    This  result  is  not  surprising  from  the 
manner  in  which  l/L  enters  U£Q  Eh 33.   A  brief  summary  of  some  of 
the  other  statistics  obtained  are  listed: 

UUo  =18707  ens       Z<,  =  0.\0U2cn     L---H?J73.bcn 
«T    =3  7847crts      So  =  0.I0238cm    ft  =  HO^fiZ.Un 

r3l.z--o.oo%\       r»i--o.ftl5'3 
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The  most  surprising  fact  regarding  the  statistics,  and 
equation  (34)  is  the  complete  unimportance  of  the  variable  Z     in 
describing  U2Q«    At  first  glance,  equation  (33)  would  suggest  a 
much  stronger  dependence,  since  -lnZ0  is  one  of  the  terms  of  the 
right  side,  and  this  should  indicate  that  U20  decreases  with 
increasing  Z0.    This  is  a  result  presented  by  many  writers,  who 
emphasize  the  tendency  for  u  to  decrease  with  increasing  surface 
Rossby  number  Vg/fZ0.     However,  from  equation  (33)  it  can  be 
seen  that 

where  (u*/k)n  is  the  neutral  value  of  (u  A).    That  part  of  (u*/k) 
which  depends  upon  stability  is  primarily  incorporated  in  the  ratio 


so  that  the  neutral  part  of  (u  A)  is  essentially 


m/c 


IH 


The  final  point  to  be  emphasized  here  then  is  as  follows: 
Within  equation  (36)  the  term  -lnZQ  tends  to  reduce  U2q-     However, 
carefully  measured  results  of  Hay  (1955)  £5]  indicate  that 


(\f\~-  I3.IJ 
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* 
so  that  u     may  be  expected  to  increase  in  the  boundary  layer  as  ZQ 

increases.     Obviously  these  two  effects  oppose  one  another,  with 

a  wind-increasing  tendency  from  (u*/k)  and  a  wind-decreasing 

tendency  arising  from  -lnZ0  within  the  bracket  of  (36).    This 

oddity  resulting  from  opposite  effects  produced  by  an  increase  in 

Z0  appears  to  explain  the  lack  of  correlation  of  U20  with  Z0. 

No  statistical  analysis  involving  u  /k  with  ZQ  has  been  made 

but   their  values      are  listed  in  the  Appendix  and  may  serve  as  the 

basis  of  future  research  on  the  subject. 
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APPENDIX  I 


Table  of  selected  micrometeorological  parameters  in  the  surface  layer. 


AT(°C)   Ri'(10"8) 


Ri 


u*/k 


1081.5 

-3.88 

-1.15 

-25197.329 

-.0057 

112.510 

.098703 

412.0 

-4.44 

-9.02 

-4900.552 

-.0327 

56.065 

.250295 

309.0 

-1.66 

-6.00 

-6399.332 

.0243 

135.457 

. 34632 

1030.0 

-3.88 

-1.26 

-23040.770 

-.0063 

107.509 

.060164 

1339.0 

-2.22 

-0,43 

-64212.942 

-.0022 

171.079 

.709676 

1081.5 

-2.77 

-0.82 

-34474.195 

-.0041 

111.440 

.01238 

1442.0 

-2.77 

-0.46 

-59822.933 

-.0024 

146.968 

.145697 

515.0 

-1.11 

-1.44 

-20410.493 

-.0071 

54.025 

.126881 

1442.0 

-6.11 

-1.01 

-28242.996 

-.0050 

149.440 

.06137 

1133.0 

-2.77 

-0.75 

-37639.751 

-.0038 

116.493 

.042618 

618.0 

-2.22 

-1.98 

-15405.752 

-.0095 

65.760 

.095000 

1133.0 

-1.11 

-0.29 

-92272.089 

-.0015 

114.868 

.111753 

2060.0 

-2.22 

0.18 

-150915.860 

- . 0009 

1 

208.067 

.158484 

875.5 

-4.44 

-1.98 

-15426.860 

-.0095 

93.154 

.089827 

772.5 

-  .55 

-0.32 

-85614.909 

-.0017 

78.378 

.098703 

618.0 

-1.11 

-0.99 

-28764.428 

- . 0049 

64.009 

.035126 

721.0 

-3.88 

-2.55 

-12423.038 

-.0119 

77.801 

.208904 

875.5 

=  1.66 

-0.74 

-37875.237 

-.0038 

55.825 

.174716 

1133.0 

-2.79 

-0.74 

-37980.775 

-.0038 

116.468 

.098703 

1030.0 

-2.22 

-0.71 

-39240.355 

-.0036 

105.800 

.098703 

721.0 

-3.33 

-2.18 

-14163.219 

-.0103 

98.529 

.738132 

1133.0 

-0.56 

-0.15 

-184229.295 

-.0008 

104.304 

.062660 

1133.0 

-3.89 

-1.02 

-29043.613 

-.0051 

117.4444 

.098703 

1287.5 

-2.78 

-0.56 

-49075.979 

-.0029 

147.914 

.315771 

1081.5 

-1.11 

-0.32 

-85030.331 

-.0017 

109.737 

.018613 

1236.0 

-1.11 

-0.25 

-110239.246 

-.0013 

125.113 

.057789 

1545.0 

-6.11 

-0.86 

-32759.215 

- . 0044 

148.858 

.087884 

772.5 

-1.11 

-0.63 

-44052.566 

-.0032 

79.153 

.098703 

618.0 

-1.67 

-1.48 

-19944.649 

-.0073 

37.485 

.023398 

721.0 

-3.89 

-2.54 

-12462.568 

-.0118 

65.775 

.103646 
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